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Abstract 

The purpose of this paper is to show that if the spectral radius of the 
matrix of the linearized system of a nonlinear discrete dynamical system is 
less then one then the characterization theorem from 2 3 holds. 



1 Introduction 

We consider the system of difference equations 

x k+1 = g(x k ) k = 0,1,2... (1) 

where g : £1 — > £1 is an analytic function defined on a domain SI included in M n . 
A point x° € £1 is a fixed point for the system JJ) if x° satisfies 

x°=g(x°) (2) 

The fixed point a: of (JIJ is "stable" provided that given any ball B(x°,e) = 
{x G Q/\\x -x°\\ < e}, there is a ball B(x°, S) = {x G Q/\\x - x°\\ < 6} such that 
if x e B(x°,5) then g k (x) G B(a; ,e), for k = 0,1,2,... gj. 

If in addition there is a ball B(x°, r) such that g k (x) — > x° as k — > oo for all 
x G B(x°,r) then the fixed point x° is "asymptotically stable". g]. 

The domain of attraction DA(x°) of the asymptotically stable fixed point x° is 
the set of initial states x G 51 from which the system converges to the fixed point 
itself i.e. 

DA(x°) = {xe n\g k (x) x } (3) 

It is known that x° is a fixed point for system (JIJ if and only if G K™ is a 
fixed point for the system 

Vk+i = f{Vk) fc = 0,l,2... (4) 
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where / : ft — x — ► Cl — x is the analytic function defined by 

f(y) = g(y + x°)-x° ioi- y en~x° (5) 

The fixed point a; of is asymptotically stable if and only if the fixed point 
6 W 1 of the system Q is asymptotically stable. 

The domain of attraction of x°, DA(x°) is related to the domain of attraction 
of 0, DA(0) by the equation 

DA{x°) = DA{0) + x° (6) 

For the above reason in the followings instead of the system Q we will consider 
the system (@J. 

Theoretical research shows that the DA(0) and its boundary are complicated 
sets. In most cases, they do not admit an explicit elementary representation. 
For this reason, different procedures are used for the approximation of the DA(0) 
with domains having a simpler shape. For example, in the case of the theorem 4.20 
pg. 170 0] the domain which approximates the DA(0) is defined by a Lyapunov 
function V built with the matrix d$f of the linearized system in 0. In we 
presented a technique for the construction of a Lyapunov function V in the case 
when the matrix <9 / is a contraction, i.e. ||(?o/|| < 1. In this paper, we extend this 
result for the more general case when r(9 /) < 1 (where r denotes the spectral 
radius). The Lyapunov function V is built using the whole nonlinear system, not 
only the matrix dof. V is defined on the whole DA(0), and more, the DA(0) is 
the natural domain of analyticity of V. 



2 The result 

Let be / : O — > O an analytic function defined on a domain fl C R™ containing 
the origin £ R". 

Theorem 1. // the function f satisfies the following conditions: 

/(0) - (7) 



r(dof) < 1 (8) 

then is an asymptotically stable fixed point. DA(0) is an open subset of O and 
coincides with the natural domain of analyticity of the unique solution V of the 
iterative first order functional equation 



V(f(x))-V(x) = -\\xf 

V(0) = w 

The function V is positive on DA(0) and V(x) — > +oo 7 for any x £ FrDA(0) 
(FrDA(0) denotes the boundary of DA(0) ). 
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Proof. Let be A = dof and r = r(dof). 

The fact is an asymptotically stable fixed point is proved in The fact 
that -DA(O) is an open subset of SI follows from the continuity of the function f k , 
for any fc G N. 

oo 

We will prove that the series ||/ fc (x)|| 2 is convergent for any x G DA(0). 

k=0 

We can decompose the function / as follows: 

f(x) = Ax + h{x) (10) 

where h : SI — > SI is an analytical function which satisfies h(0) = 0. 
As r < 1, there exists c > such that 

|| A k x\\ < cr k \\x\\ Va:efi,fceN (11) 

Let be c = max{l,c} and e — fe? > 0. As /i is continuous and ft.(0) = 0, it 
follows that there exists 6 > such that 

||/i(x)|| <e||x|| VseB(0,i). (12) 

Let be x G DA(0) and Xk = f k (x). This provides the existence of k x G N such 
that Xfe € -6(0,5), for any k > k x . Let be yk = ^fc+fc^- This sequence satisfies 
Uk G 5(0, (5) for any fc G N and it also satisfies 

The formula of variation of constants gives: 

fe-i 

y fe = ^ fe 2/o + ]T ^-^hiyi) Vfc G N* (13) 

i=! 

Relations Ijllfl and (|13fl provide 

fc-i 

llMll ^c^llyoll+^cr^- 1 ^^)!! Vfc g N* (14) 
t=o 

and using l|12f> and that c < c, the following inequality follows: 

fc-i 

\\y k \\<cr k \\y \\+Y l cr k - i - 1 e\\y i \\ Vfc G N* (15) 

i=0 

Relation l|15(l can be written as 

k-l 

r- k \\ V k\\ < c\\Vo\\+J2 ^-^(r-ij/ill) VfcGN* (16) 

i=0 

Gronwall's inequality for the discrete case provides 

fe-i 

r- fc ||y fc || <c|| 2/0 |in(l + ^- 1 s) VfcGN* (17) 

i=0 
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thus 

llVfcll < r fe c||y ||(l + cr^ef = c\\y \\(r + ce) k = c\\y \\(^±±)» Vfc e N* (18) 
Denoting £ii = a < 1, the relation 1|18[1 gives 

\\x k \\ < c\\x k Ja k Vk>k x (19) 
Thus, for any x G DA(0) there exists fc x > such that 

\\f k (x)\\<c\\f k *(x)\\a k Vk>k x (20) 

oo 

which assures that the series ^2 \\f k (x)\\ 2 is convergent for any x G DA(0). 

k=0 

Let be V = V(x) the function defined by 

oo 

V(x) = J2\\f k (x)\\ 2 VxeDA(O) (21) 

fc=0 

The above function defined on DA(0) is analytical, positive and satisfies ©. In 
order to show that the function V defined by l|21|l is the unique function which sat- 
isfies (@ we consider V' = V'(x) satisfying 10 and we denote by V" the difference 
V" = V - V. It is easy to see that V"{f{x)) - V"(x) = 0, for any x G DA(0). 
Therefore, we have V (x) = V"(f k (x)) for any x G DA(0) and any k G N. It 
follows that V"{x) = lim V"(f k {x)) = for any x G DA(0). In other words, 

k — >oo 

V(x) — V'(x), for any x G DA(0), so V defined by l|21|l is the unique function 
which satisfies ©• 

In order to show that V(x) qq f or an y x° G FrDA(0) we consider x° G 

FrDA(0) and r > such that ||/ fe (a; )[| > r, for any k = 0,1,2.... For an 
arbitrary positive number AT > we consider the first natural number k\ which 
satisfies k\ > + 1. Let be r 1 > such that ||/ fc (x)|| > 4- for any 

and x G B(x°,ri). For any a; G B(x°,n) n DA(0) we have £ ||/ fc (a;)|| 2 > N. 

k=a 

x^x" 

Therefore, V(x) — > oo. □ 



3 Numerical examples 

In this section, some examples of discrete systems are given, for which the domain 
of attraction of the zero steady state can be estimated using the method described 
in El 



Example 1 

The following system of difference equations is considered: 



Xn+i x n y n -f" y n 

Vn+l = Vn 
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It is clear that (0,0) is an asymptotically stable fixed point for the system 1)22(1 . 
It can be proved theoretically that the domain of attraction of (0, 0) is DA(0) = 
R x (-1, 1). 

The matrix of the linearized system has the norm ||<9o/|| = 1, thus the theorem from 
does not apply. On the other hand, r(dof) = 0, thus the characterization 
theorem for DA(0) given above applies. Thus, the domain of attraction can be 
estimated using the numerical method given in The first step of the method 
gives the whole domain of attraction: 
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Figure 1: The estimate of DA(0) after 1 step for system l(22|l 
Example 2 

The following system is considered: 



x n +i = 4x„ 4 

Vn+1 — ^nUn 



Vn 



(23) 



This system of difference equations has the asymptotically stable fixed point (0, 0). 
It is once again the case when ||3o/|| — 1 an d K^o/) = 0- After one step, we can 
obtain the following estimate of the DA(0): 
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Figure 2: The estimate of DA(0) after 1 step for system <|2l?|l 
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